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ABSTRACT 

Disk galaxies at high redshift [z ^ 2) are characterized by high fractions of cold gas, strong tur- 
bulence, and giant star-forming clumps. Major mergers of disk galaxies at high redshift should 
then generally involve such turbulent clumpy disks. Merger simulations, however, model the ISM 
as a stable, homogeneous, and thermally pressurized medium. We present the first merger simula- 
tions with high fractions of cold, turbulent, and clumpy gas. We discuss the major new features of 
these models compared to models where the gas is artificially stabilized and warmed. Gas turbu- 
lence, which is already strong in high-redshift disks, is further enhanced in mergers. Some phases 
are dispersion-dominated, with most of the gas kinetic energy in the form of velocity dispersion and 
very chaotic velocity fields, unlike merger models using a thermally stabilized gas. These mergers can 
reach very high star formation rates, and have multi-component gas spectra consistent with SubMil- 
limeter Galaxies. Major mergers with high fractions of cold turbulent gas are also characterized by 
highly dissipative gas collapse to the center of mass, with the stellar component following in a global 
contraction. The final galaxies are early-type with relatively small radii and high Sersic indices, like 
high-redshift compact spheroids. The mass fraction in a disk component that survives or re-forms 
after a merger is severely reduced compared to models with stabilized gas, and the formation of a 
massive disk component would require significant accretion of external baryons afterwards. Merg- 
ers thus appear to destroy extended disks even when the gas fraction is high, and this lends further 
support to smooth infall as the main formation mechanism for massive disk galaxies. 

Subject headings: galaxies: formation — galaxies: interactions — galaxies: high-redshift — galaxies: 
elliptical and lenticular, cD — galaxies: structure 



1. INTRODUCTION 

In the A-CDM cosmological model, collisions and 
mergers are an important growth mechanism for galaxies, 
and most o f them should have occurred at high redshift 
(z > 1, e.g.. lStewart et al]l2009f ). Major mergers, where 
the ratio of the baryonic masses of the involved galax- 
ies is close to one, are in particular important for the 
build-up of massive galaxies (Hopkins ct al. 2010). Even 
if recent models and high-redshift observations suggest 
that the growth of galaxies near i* is largely through 
relatively smooth and cold accretion (jDekel et al.l 120091 : 



Bournaud fc Elmcgrccn"2009t IConrov fc Wechsleri 120091: 
Genel et al.ii2010t) , about one-third of a galaxy's baryons 
are s till expected to be provided by signifi cant mergers 
(e.g., iDekel et al.|[20Ql iBrooks et aLll2009l ). 

The effect of the merger of two disks galaxies has 
been largely explored in the low-redshift context, where 
galactic disks are mostly made up of stars with relatively 
modest fractions of cold gas. In such conditions, it 
is known that a major merger of two spiral galaxies 
typically produces a remnant re sembling an ellipti- 
cal or early-type galaxy (ETC) ()Barnes fc Hernqtiis^ 



19961: iMihos fc Hernquisti 119961: iNaab fc Burkerl 
200a iBournaud et all 2004al: iNaab et all 120061 120071: 



Bournaud et al.ll2007bl: iHoffman et al.ll2010[ ). There are 



still some differences bet ween the predictions of merger 
models and observations (jBurkert et al.l 12008^ but there 
is, overall, an ample consensus that hierarchical merging 
is the main formation channel for present-day ETGs. 

At high redshift, the role of major mergers in the 
formation of particular galaxies and in the growth 
of galaxies in general remains much more debated. 
First, determinations of the merger rate from high- 
redshift surveys are usually based on irregular struc- 
tures or disturbed kin e matics Ce.g.. iConselice et alj 
200l ILotz et alJ [20061 [20081 lOverzier et alJ 120101: 



Jogee et alJ 120091: iLopez-Saniuan et alJ 120091: iPuechI 
20101 ). However, primordial galaxy disks can have 



irregular morphologies and disturbed kinematics even 
when they are mos tly isolated, not undergo- 
ing strong interactions ([Elmegreen et al.l 120071 120091: 
Genzel et alll2008t IBournaud et al.ll2007al: iGenzel et al.l 



20081: IBournaud et all 120081: iDekel et al.l |2009[) . so 

the actual merger rate remains unknown. Also 
high-redshift disk galaxies have ty pical morph o logies 
that differ from nearby spi r als (jC ow ie et al.l 11996 



Ehnegreen fc ElmegreenI 120051 : lElmegreen et al 



2007 



Genzel et al.ll2008D . and their mergers could also have dif 



ferent morphologies and kinematics compared to nearby 
mergers, potentially making the identification of mergers 
more ambiguous: for instance, asymetries and dumpi- 
ness in_0£tka]_jmagingare often attributed to merg- 
ers (jConselice et al.l 120031: ILotz et al.l [2008D . but high- 



redshift disk galaxies seem to frequently form massive 
clumps just b y internal instabilities that do not require 
intera ctions (jElmegreen et al.l 120071 120091 : IGenzel et all 
120081 ). Even in cosmological simulations, the galaxy 
merger rat e is continuously revised (e.g., iGenel et al.l 
120081 [20T0I and references therein). In particular, merg- 
ers of dark halos do not necessarily lead to the merging 
of their central galaxies, and the baryonic mass ratio of 
galaxy mergers can significantly differ fro m the dark mat - 
ter mass ratio of their host haloes (Hop kins et al.l [201 Of ). 

Second, the outcome of major mergers is usually con- 
sidered to be early- type galaxies (ETGs). However, high- 
redshift ETGs are often unexpectedly compact compared 
to ne arby ellipticals and r uerger simulation s (Daddi ct ahl 
2005"; 'Truiillo et al.l|2006|; Ivan Dokkum et al.ll2009l l200a 
Krick ct al. 2009). The role of mergers in forming 
these ETGs has been questioned: while their evolu- 
tion into modern, more extended elliptical galaxies can 
be explained thro ugh continuous hierarchical merging 
(jNaab et al.ll2009[ ). the emergence of compact ellipticals 
at high redshift remains more mysterious. 

Third, the most general issue is the ability of disk- 
dominated galaxies to survive violent mergers. While 
it is known that mergers of disk galaxies can produce 
ETGs with no or little residual disk component (ref- 
erences above), whether all mergers even with very 
high gas fractions destroy disk galaxies and transform 
them into ETGs, or whether disk-dominated galax- 
ies could survive major mergers in high-redshift con- 
ditions, remain more uncertain. The frequency of vi- 
olent mergers could indeed call i nto question the sur - 
vival of disk-dominated galaxies (jWeinzirl et al.|[2009D . 
Some models of mergers, however, have led to the 



proposal that high-redshift merger remnants could be 
dominated by large rotati ng disk components rather 
than bulges and sphe roids ([Springel fc Hernquisti 120051 : 
[Robertson et al.l|2006[) . There have been claims that such 
disky merger remnants cannot account for the ob served 
properties of high-redshift disk gala xies (Shapiro et al.l 
120081: IBournaud fc ElmegreenI 120091) but this remains 
actively debated in simul ations ([Robertson fc Bullockl 
I2OO8D and in observations ([Hammer et al.ll2009f ). 

Thus, at least three major unknowns remain: the 
frequency and observational signatures of high-redshift 
mergers, their role in forming the high-redshift popula- 
tions of ETGs, and the ability of massive disk-dominated 
galaxies to survive these events. 

The dynamics and outcome of mergers have never been 

modeled for realistic , high-redshift galaxie s. Simula- 

tions by [Cox (200 41). iSpringel fc Hernq"^ (|2005[ ). and 
[Robertson et al.l (|2006f ) considered high gas fractions and 
concluded that disks could survive or re-form after a 
merger. However, these simulations had high thermal 
pressure support in the ISM (temperature>10''K), which 
forces the gas to be smooth and stable in the isolated 
galaxies and in the colliding pairs as well. 

Star-forming galaxies at high redshift (z = 1 — 
5) are very different from such models. They 
are generally gassy, clumpy, and turbulent, with 
^ 50% of their baryons in cold m olecular gas 
([Daddi et al.l l2010al: iTacconi et all [20 10D, giant star- 
forming clu r ups o f 10^ — 10® Mq ([Elmegreen et al.l 
[2004 [20071 I2009D . and velocity dispersions of sev- 
eral tens of km s~^ (iForster Schreiber et al.l 120091: 
Shapiro et aH l2008t fGenzel et al.l 120081 : lEpinat et al.l 
20091 iQuevrel et al.l |2009[ ). Associated star formation 



rates (SFRs) are around 100 Mq yr'^. While 
most observations of high-redshift star forming galaxies 
concerned only very massive galaxies (above the crit- 
ical mass L*), recent studies of z > 1 lensed galax- 
ies find that gravitational instabilities and turbulent 
clum ps could also be dominant in lower-m ass galax- 
ies ( Jones et al.l l2010t iSwinbank eFall |2010[ ). The 
morphology and dynamics of the disks suggest that 
the observed clumps fo r m by gravitational instabilities 
JElmegreen et al.' '2004 1: lElmegree n fc Elmegreen' '20051 
Shapiro ct al. 2008; Bournaudet al. 2008; Daddi ct aL| 
2010ai) . Simulations of this process require cosmological 



codes that allow for highly supersonic tur bulent speeds 
([Agertz et al.l [20091: ICeverino et al.l 120101 ). Even low- 
redshift spiral galaxies have an unstable, cloudy ISM 
supported mostly by s upersonic turbu lence rather than 
thermal pressure (e.g., iBurkertI I2006D . Recent simula- 
tions showed that a spatial resolution of 10% of the scale 
height with full tracking of turbulent motions is required 
to reproduce the three-dimensio nal inertial range of tu r- 
bulence inside a modern galaxy ([Bournaud et al.ir2010[ ). 

A typical wet merger at high redshift should involve 
the interaction of such clumpy, turbulent disks. The 
outcome could differ from previously modeled merg- 
ers using thermally-stabilized gas. Recent studies have 
shown that ISM turbulence and dumpiness have a sub- 
stantial effect in mergers of present-day spirals with 
just a few percent of gas, with significant differences 
both in the star-formation history during the interaction 
([Tevssier et al.. ,2010; .Saitoh et al.. .2009.) and the prop- 
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erties of the relaxed p ost- merger ETG (jBournaud et alj 
120081 iBois et alJlMol) . The effect could presumably be 
more dramatic in high-redshift mergers involving high 
fractions of cold gas. 

In this paper, we present hydrodynamical AMR sim- 
ulations of major mergers between galaxies with high 
fractions of cold, turbulent and clumpy gas, like typical 
z ~ 2 star-forming galaxies. Gas cooling below 10^ K is 
allowed, supersonic ISM turbulence is captured, and the 
main star-forming complexes are directly resolved. We 
compare with models using artificially stabilized disks, 
as in traditional high-redshift studies. 

We find that mergers of high-redshift disks can un- 
dergo a very clumpy and chaotic phase, during which 
the kinematics is dispersion-dominated even for the gas 
component. The gas kinematics appears consistent with 
the observed properties of some "dispersion-dominated 
galaxies" as well as the spectral properties of SubMil- 
limeter Galaxies (SMGs) . A large part of the total mass 
is thus supported by supersonic turbulence, which is a 
dissipative support, hence the merging system undergoes 
a rapid dissipational collapse into a compact ETG. The 
mass fraction in a surviving or re-formed disk component 
is much lower for supersonically turbulent high-redshift 
mergers than for models with artificially warmed gas. 
This does not mean that mergers of gas-rich galaxies can 
never re-form a disk: a large enough far-outer reservoir 
of gas could recollect into a new disk, or continuous in- 
fall of fresh gas could build a new disk. Some orbital 
parameters not modeled here might also be m ore favor- 
able to disk survival (e.g. iHopkins et ani2009l ). Still, our 
main conclusion is that mergers between gas-rich turbu- 
lent and clumpy galaxies are qualitatively different than 
mergers between gas-rich smooth galaxies with artificial 
pressure support, and can better explain resolved obser- 
vations of high-redshift systems. 

2. SIMULATIONS 

This paper studies idealized models of disk galaxies. 
The initial conditions are not cosmological, and there is 
no cosmologically-motivated boundary condition model- 
ing mass infall or extended gas reservoirs. Instead, the 
initial conditions are designed to be representative of gas- 
rich major mergers of redshift z ~ 2 galaxies, based on 
their main observed characteristics. We model collisions 
and mergers between two such galaxies and study the 
effect of this merger, by comparison to the modeled evo- 
lution of a single isolated galaxy. 

The sim ulations were p erformed with the AMR code 
RAMSES (|TevssieH[200l . The technique is fully de- 
scribed elsewhere and shown to model realistic inter- 
stellar gas in low-redshift mergers and is olated disks 
(|Tevssier et al.l[2010l : iBournaud et al.ll201 0^. The main 
parameters adopted for the present simulations are in 
Table 1. 

The box size for the simulation is 200 kpc. The coarsest 
level of the AMR grid is / = 8, which corresponds to a 
256"^ Cartesian grid with a cell size of 781 pc. A cell 
that contains a gas mass larger than TOj-os = 8 x lO"* Mq, 
or a number of particles larger than 15, is refined, until 
the maximal level / = 11 is reached. At that point, the 
cell size corresponds to 97 pc. Stars and dark matter are 
described with 2 x 10^ coUisionless particles each in the 
initial galaxies. 



Our clumpy disk simula tions were performed with a 
barotropic cooling model (jBournaud et al.l [20101 ) . nat- 
urally producing a cloudy and turbulent ISM, with a 
temperature floor around 10'^ K. We refer to these as 
the cooling models. For comparison, the same gas-rich 
mergers were modeled with a thermally pressurized and 
Toomre-stable ISM, using an adiabatic Equation of State 
(EoS) with an exponent of 7 = 5/3 for densities above 
1 cm^^, and r=10'' K for lower densities. This EoS 
maintains a Toomre parameter Q ~ 1.5 — 2 in the initial 
disks and stabilizes the gas against axisymmetric per- 
turbations. These will be called the stabilized models. 
In all density-dependent pressure floor ensures 

that the Jeans length is resolved by at least 4 cells to 
avoid artificial fragmentation, as initially proposed by 
Machacck et al. (2003) (see Teyssicr et al. 2010 for the 
implementation details) . 

In the cooling models, star formation is assumed to 
proceed above a density threshold of 300 cm~^ with an 
efficiency of 7%, i.e. 7% of the gas in a given cell forms 
stars per local free-fall time. This gives a star formation 
rate (SFR) of 120 Mq yr"^ in the isolated disk model, 
realistic for such gas-rich z ^ 2 systems. In the stabi- 
lized models, the star formation threshold is 3 cm~^ and 
the efficiency is 4%, giving the same SFR in the isolated 
disks: we compare the two sets of models with the same 
SFR in pre-merger disks. The resulting star formation 
history in the merger models with cooling/stabilized EoS 
is relatively similar (see Section 1X4)) . 

All sim ulations use the kinetic f eedback model de- 
scribed in [Dubois fc Tevssieri ()2008[ ). with 20% of the 
energy from each supernovae re-injected in the form of 
an expanding bubble of initial radius 100 pc (assuming 
that the rest of the supernovae energy was radiated away 
before propagating up to scales of 100 pc). Some mod- 
els have increased feedback with 100% of the supernova 
energy injected into the ISM (models labeled "F") and 
some have a supernova feedback efficiency set to 500% 
(models labeled "F5", discussed only in Section 3.4). 

Our simulations start with disk galaxies, each having 
a stellar mass of 4 x 10^° Mq, an initial gas fraction of 
70%, a disk scale length of 5 kpc with a truncation radius 
of 12 kpc, and an initial scale height of 800 pc. An initial 
bulge containing 15% of the stellar mass is assumed, with 
a Hernquist profile and a scale length of 500 pc. The dark 
matter halo has a Burkert profile with a core radius of 
8 kpc, and a mass fraction (dark/total) inside the disk 
radius of 30%. The circular velocity in the outer disk is 
245 km s~i. 

In the cooling models, the initial disks spontaneously 
become clumpy and turbulent, with V/a ~ 5 (for rota- 
tion speed V and turbulent speed a). When the mergers 
occur, at the first pericenter, gas fractions are around 
50% after some early gas consumption (s ee Table 2): 
this seems typical for disk galaxies at z ~ 2 (jPaddi et al.l 
l2010at iTacconi et al.ll2010l ) , and hence expected for a wet 
merger at z ~ 2. 

Our merger models start with relatively massive disk 
gal axies, so as to c o mpar e with existing models, e.g. 
in [Robertson et al.l ([2006[ ). There are galaxies like 
th is in z ^ 2 samples (for instance in the sample 
of [Forster Schreiber et all [20091 ). Also, our pre-merger 
disks are not extreme in terms of gas dumpiness and 
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TABLE 1 

Parameter used for the simulations. 



Model 


Orbit 


EoS 


Comment 


Oi 


1 
i 


cooling 


merger of gas-rich clumpy disks 


GIF 


1 


cooling 


CI with stronger feedback 


C1F5 


1 


cooling 


CI with 500% feedback 


SI 


1 


stabilized 


gas-rich, smooth, stabilized disks 


SIF 


1 


stabilized 


SI with stronger feedback 


S1F5 


1 


stabilized 


SI with 500% feedback 


02 


2 


cooling 


merger of gas-rich clumpy disks 


S2 


2 


stabilized 


gas-rich, smooth, stabilized disks 


C3 


3 


cooling 


merger of gas-rich clumpy disks 


S3 


3 


stabilized 


gas-rich, smooth, stabilized disks 


I 


isolated 


cooling 


control run: high-redshift disk 


LM-C2 


2 


cooling 


lower-mass gas-rich clumpy disks 


LM-C2F 


2 


cooling 


C2 with stronger feedback 


LM-S2 


2 


stabilized 


lower-mass, smooth stabilized disks 



TABLE 2 

Properties of the final relaxed systems: 

Half-mass radius R\/2 a^nd Sersic index n (best fit between 0.3 and 3 R\/2)\ g^s fraction remaining at the pericenter and 400 Myr later; 
mass-weigthed average of the ID local gas turbulent speed < (Jtuih.iu > measured 140 Myr after each pericenter; baryonic mass fraction 
of the baryons in a rotating disk /diski measured using a kinematic disk-f spheroid decomposition as in Martig & Bo urnauSI I I2010I) . 





Gas consumption 


ISM turbulence 


Merger remnant properties 


Model 


/gas, peri 


/gas,i-|-400Myr 


< o"turb,iD > (kms^i) 


Sersic n 


Rl/2 


/disk 


CI 


45% 


17% 


175 


4.4 


2.8 


12% 


GIF 


54% 


22% 


186 


4.2 


3.0 


16% 


C1F5 


58% 


21% 


204 


3.9 


3.1 


21% 


SI 


48% 


20% 


56 


2.9 


6.1 


40% 


SIF 


52% 


22% 


59 


2.5 


6.3 


48% 


S1F5 


56% 


23% 


65 


2.6 


6.2 


53% 


C2 


53% 


13% 


153 


4.2 


2.9 


14% 


S2 


49% 


18% 


68 


3.2 


5.8 


37% 


C3 


53% 


19% 


164 


4.9 


2.6 


11% 


S3 


47% 


17% 


72 


3.3 


5.6 


37% 


I 






59 


1.7 


5.4 


66% 


LM-C2 


52% 


19% 


89 


3.7 


1.2 


10% 


LM-C2F 


58% 


21% 


96 


3.9 


1.3 


13% 


LM-S2 


57% 


22% 


42 


2.7 


2.0 


35% 



Gas - CI 




T=-155Myr 



Gas - S1 





t=+55Myr 



30kp<; 




V 



1=+215Myr 



t=+360Myr 




Fig. 1. — Snapshots of the gas mass distribution for models CI (cooling) and SI (stabilized ISM) at similar instants and under similar 
projections, t = is the first pericenter passage and all maps are in log scale. 
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Fig. 2. — Snapshots of the stellar mass surface density distribution for models CI (cooling) and SI (stabilized ISM) at similar instants 
and under similar projections, i = is the first pericenter passage and all maps are in log scale. 
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Fig. 3. — Three orthogonal projections of model CI showing the 
surface density, velocity and dispersion maps (mass-weighted val- 
ues along the line-of-sight, for gas denser than 10 cm""^). A res- 
olution of 1 kpc was assumed (FWHM of a gaussian beam). All 
snapshots are 44 X 44 kpc. Note that Ho observations would be 
mostly sensitive to emission from the dense clumps. 

turbulenc e: many z ^ 2 disks have ev en lower V/a ra- 
tios (e.g.. iForster Schreiber et aI1l2009D . 

Three interaction orbits were used. They are all pro- 
grade for one galaxy and retrograde for the other in or- 
der to consider a typical total angular momentum rather 
than extreme cases with aligned spinfQ. Orbit 1 has an 
impact parameter of 12 kpc and is parabolic. One disk 
has an initial inclination of 30° from the orbital plane, 
and the other has an initial inclination of 50°. For or- 
bit 2, these parameters are: 35 kpc, hyperbolic (total 
energy of the galaxy pair is 0.3 times its initial kinetic 
energy), and 30° and 50° degrees, respectively. For or- 
bit 3, they are 25 kpc, hyperbolic (total energy = 0.2 x 
the initial kinetic energy), 40°, and 65° degrees, respec- 
tively. 

Three simulations were also run for a merger of lower- 
mass galaxies. The masses of all initial components (stel- 
lar disk, gaseous disk, bulge, halo) where reduced by a 
factor of 8 compared to the initial higher-mass models, 
and all sizes and distances were reduced by a factor of 
^/8, hence keeping surface densities about constant. The 
pre-merger galaxies then had circular velocities of about 
150 km s~^ in their outer disks, with a F/tr ~ 3.5 ra- 
tio - somewhat lower than in the more massive disks 
models, which is in agreeme nt with observed trends 
(|Forster Schreiber et alj I2009D . These merger models 
were performed using Orbit 2 with our cooling and sta- 
bilized ISM models, and are respectively denoted LM-C2 

^ Aligned spins would not necessarily result in higher angular 
momentum or more disk component in the final result, because of 
tidal removal of the angular momentum (Robertson et al. 2006 ) 
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Fig. 4. — Same as Fig. 2, with comparison of on-going mergers in 
the cooling and stabilized ISM models at the same instant and un- 
der the same projection, shown here at full resolution. An isolated 
clumpy turbulent disk is also shown after a similar evolutionary 
time as the mergers. 

and LM-S2, and LM-C2F with increased feedback. 

3. RESULTS 

3.1. On- going mergers 

The time evolution of mergers CI (cooling) and SI 
(stabilized ISM) is shown in Figures [1] and [5] The sta- 
bilized model has relatively homogeneous tidal tails or- 
biting around smooth gas disks. The gas in the cooling 
model is dominated by numerous star-forming clumps 
that are somewhat more massive than the clumps were 
in the pre-merger galaxies. Visual inspection shows that 
some of the pre-existing clumps survived and some have 
disrupted and re-formed. Numerous filaments form in 
many directions, instead of a main pair of continuous 
tidal tails. Other projections (Fig. [S]) show that the gas 
distribution during the merger forms a very irregular 
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radius (kpc) radius (l<pc) radius (kpc) 

Fig. 5. — Radial profiles of the stellar surface densities in final merged and relaxed systems, and the isolated disk model after the same 
evolutionary time. Several projections were averaged for each model. Ellipse-fitting was performed with the ellipse task in IRAF. 



spheroid in three dimensions; there are no disks. The 
model with stronger feedback (GIF in Fig. 3) is quaU- 
tatively similar: self-gravity is already sufficient to drive 
strong ISM turbulence in massive galaxies at high red- 
shift (Dckcl et al. 2009; Elmcgrcen & Burkcrt 2010). 

Gas velocity fields and line-of-sight dispersion maps are 
shown in Figures [3] and S) The gas velocity fields of our 
high-redshift cooling models are very chaotic, and often 
lack extended rotating components, such as large disks or 
long tidal tails with monotonic velocity gradients. The 
gas velocity dispersions are high, especially near dense 
clumps. 

Interactions at low redshift substantially increase the 
gas velocity dispersion, from typically 10 km s~^ in non- 
interacting spirals to 3 0-40 km s~^ in m ajor mergers (see 
Elme green et al.l [19951 for observations, iBournaud et al.l 
20081 for simulations). We find a similar relative increase 
here, but starting with disks that are already quite tur- 
bulent before the merger. This results in systems where 
the gas component is dispersion-dominated during the 
merger, with V j a ratiofQ of around 2 or even 1 for some 
projections and times. The gas dispersions are high in 
projection on the line-of-sight, and also in direct mea- 
surements of the local three-dimensional motions: we 
measured the turbulent velocity dispersion in 1 kpc'^ 
cubes in model Gl at t — 140 Myr; the mass- weighted 
average value was < crturb.iD >— 175 km s~^ (higher 
values on projected maps such as Figures 3 and 4 result 
from linc-of-sight projection effects, not to a physical, lo- 
cal turbulent motion). Similar measurements are given 
for all models in Table 2, at the same instant which is 
about the peak of turbulent speed in our merger models 
with cooling. These measures show that the turbulent 
velocities become much higher during the interactions 
and mergers than in the pre-merger clumpy disks, and 
that this property is not reproduced in the simulations 
using a stabilized gas model. The values of line-of-sight 
velocity dispersions in Figures 3 and 4 will be further 
compared to high-redshift observations in Section 4.1. 

Mergers with the stabilized ISM model show velocity 
dispersions that are lower by a factor ~ 3 for all times. 
Instead of chaotic velocity fields, they have relatively 
smooth and extended rotating gas disks, surrounded by 
long tidal tails that co-rotate with large-scale velocity 
gradients (see examples in Fig. 3). 

Low-redshift mergers have different signatures, since 
they typically exhibit long tidal tails with large and 

^ Throughout the paper, a refers to a one-dimensional disper- 
sion. A system with V ja < 2 is dispersion-dominated with > 60% 
of its kinetic energy support is in the form of dispersions. 



smoot h velocity gradients along these tails (see exam - 
ples in IBournaud et al.ll2004bl : iChilingarian et al.ll2010D . 

Low-redshift merger models with warm stable gas ver- 
sus colder and unstable gas in cooling models have been 
compared in iTeyssicr ct al. (2010). While they display 
substantial differences in their star formation properties, 
both models are dominated by long tidal tails and have 
similar half-mass sizes after the merger. The tails have 
star-forming clumps in the cooling model, but the clumps 
do not dominate the morphology and kinematics as they 
do at higher redshift. 

In this sense, stabilized ISM models are a more ac- 
ceptable recipe for the large-scale gas properties in low- 
redshift mergers than in high-redshift ones. This can 
be explained because low-redshift systems have weaker 
gas turbulence and are more stable than high-redshift 
systems on scales of a kpc. Low-redshift disks are only 
slightly unstable, and the sizes and masses of the clumps 
that form are much smaller than they are at high red- 
shift. Thus the effects of ISM turbulence and dumpiness 
are more striking in high-redshift mergers. Resolving 
gas turbulence, star forming instabilities and/or shocks 
change the star formation history of these lo w-redshift 
mergers (.Barnes 2004: Saitoh et al.ll2009l : iTev ssicr ct al] 
l20Ta [Cifim fc BarnesllMC but have no major impact 
on the morphology of the merger r emnants and on ly a 
limited impact on their kinematics ()Bois "eFalllmOl) . 

3.2. Merger remnants 

Supersonic gas turbulence dissipates rapidly, which can 
produce a strong contraction of a dispersion-supported 
gaseous system. Such dissipative collapse is observed in 
Figure 1 for model Gl. Some clumps are expelled by tidal 
interactions, but most of the gas clumps and inter-clump 
gas coalesce in a compact central object. Because the 
gas fraction is initially high, a large part of the final stel- 
lar content forms within the same gas clumps and follows 
their coalescence. Old stellar components contract due to 
the dissipation of the gas and young stars, which contain 
a large fraction of the total mass. Thus, the merger leads 
to the formation of a relatively compact stellar spheroid. 
This spheroid has a low half-mass radius, a high Sersic 
index, and an extended stellar halo visible as a low sur- 
face mass density component beyond radii of ^ 7 kpc 
in Figure 5. Table 2 gives half- mass radii and Sersic 
indices for all models. Figure 4 shows the correspond- 
ing stellar mass profiles. Remnants of the stabilized-ISM 
merger models are more extendecQ and have lower Ser- 

Merger remnants in the stabilized models are more extended 
in terms of half mass radii: the outer, low-mass density halo is more 
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sic indices than remnants of the coohng merger models. 
The isolated clumpy turbulent disk model has formed a 
central bulge throu gh internal e volution (as in [ Noguchil 
Il999lllmmeli et all 2004and Elm egreen et alJl2008l mod- 
els), but remains d ominated by a rotating exponential 
disk (jBournaud et al. 2007a). 

The merger remnant of a typical high-redshift merger 
is not just more compact, but also less disky. Gas disk 
re-growth is faster in the models with a stabilized ISM 
(see Fig 1 at t = 215), and thus forms a more massive 
final stellar disk. For example, the final disk component 
in model SI contains 40% of the baryons and is obvious 
up to large radius in Figure 1. In cooling model CI, 
most star formation takes place in gas clumps before they 
coalesce and get a chance to re-form a disk; the final 
disk in model CI is compact, with a disk fraction of only 
12%. Similar differences are found for all merger orbits 
(Table 2). 

The increased dissipation and resulting compactness 
in cooling models with a clumpy turbulent gas can be 
explained: clumps interact with each other gravitation- 
ally, scatter, and develop 3D motions (i.e., high disper- 
sions). The gas clumps and holes can then pass by each 
other. In stabilized ISM models, the smooth distribution 
of the gas does not allow it to pass through itself and 
it therefore rapidly settles back into a relatively planar 
disk. Gas dumpiness can reduce the ISM cross section 
for self-interaction and thereby increase the dissipation 
timescale, but the supersonic turbulence that accompa- 
nies this dumpiness is highly dissipative. For modest 
clump filling factors, the net result is a greater dissipa- 
tion rate for ISM kinetic energy. This is unlike the situa- 
tion in stabilized ISM models, where the gas is supported 
by a steady thermal pressure. Excess thermal energy 
can be radiated away, but the EoS, high temperature 
floor, and/or feedback recipe, keep the gas warm in these 
models and make the thermal support non-dissipative in 
practice. Furthermore, the massive clumps that form in 
highly turbulent gas rapidly migrate inwards through dy- 
namical friction, unlike the relatively homogeneous gas 
in stabilized models which settles into a circular orbit. 

3.3. Interpretation: increased turbulence and rapid 
dissipation in mergers 

The effect of turbulence dissipation on the size evolu- 
tion of a merger can also be estimated through simple 
calculations. A star and gas system in equilibrium satis- 
fies the Virial Theorem, and then the total energy from 
self-gravity and motion equals half of the gravitational 
potential energy, Etot — 0.5i?grav, both of which are neg- 
ative. Loss of energy by dissipation leads to a more neg- 
ative i?tot, and so a smaller, more tightly-bound equi- 
librium. For typical mass distributions in a galaxy, the 
gravitational potential energy scales with the inverse of 
the radius, _R, in which case the relative contraction is 
AR/R = -AE/E = AE/\E\. 

If the gas mass is Mgas and the ID gas velocity disper- 
sion is cr, then AE = — l.SMgasCr^ when a significant frac- 
tion of the gas kinetic energy is lost. With similar nota- 

extended in the cooling models, but would be hardly detected in 
the majority of high-redshift observations, as shown by Mancini et 
al. (2010). 



tion for the stars, the total energy is i?tot = l.SMgasO"^ -I- 
l.SMstars^^, assuming that stellar motions are of the or- 
der of the circular velocity V (even if not consisting in 
organized rotation). Thus AE/E = —X/{1 + X) where 
X = (Mgas/M,tars) X T^, with T = a/V. This ratio, 
—X/{1 + X), is also the relative contraction of the ra- 
dius after all of the initial gas motion dissipates. If fg is 
the disk gas fraction, then X — {fg/ [1 — fg]) r^. 

Turbulent dissipation sho uld occur in about a cross- 
ing time for radial motions ()McLowlll9990 . The crossing 
time for a whole galaxy is about the orbit time, typ- 
ically 100-200 Myr, so most of the non-rotational gas 
kinetic energy can be dissipated over the duration of a 
major merger. After dissipation and galaxy contraction, 
the gas receives new kinetic energy from the potential 
energy, and can then dissipate even more energy. This 
cycle of dissipation, contraction, heating and more dis- 
sipation, continues with ever decreasing orbit time until 
the radial component of the turbulent energy is gone. 

Low-rcdshift Milky Way-like disks have r ^ 3 — 
5%, and their major mergers have a val ue about 3- 
5 times higher, r ^ 20% (observatio ns: llrwini 119941: 
Elmegreen et al.l 119951 - simulations: iBournaud et al.l 
20081: iTevssier et al.ll2010t) . Together with fg of a few 
percent, the size of a merger remnant is not much smaller 
than the original galaxy size. The same would be true in 
high-redshift merger models that have subsonic gas, i.e., 
gas without much energy dissipation, even if fg is large. 

High-redshift disk galaxie s have t ~ 20% (e.g. 
iForster Schreiber et alj|2009[ ). fg ~ 0.5 (e.g., Daddi et 
al. 2010, Tacconi et al. 2010), and highly dissipative 
gas. With an increase of gas dispersions of a factor 
three during interactions (conservatively based on the 
known increased at low-redshift, but also consistent with 
high-redshift mergers spectroscopy in, e.g., Shapiro et al. 
2008), T would become of the order of ^ 60% for high- 
redshift mergers. Then AR/R ~ 0.25. This result im- 
plies that turbulent dissipation for a gas-rich system can 
cause significant radial contraction of both the gaseous 
and stellar components for every gas dissipation time, 
i.e. just 100-200 Myr. As the gas continues to dissipate 
energy, the systemic contraction of the galaxy continues, 
and the total contraction of the merger can be larger if 
some mass is expulsed from the system and carries more 
energy away. 

Another mechanism for galaxy contraction is inward 
clump migration through dynamical friction. The 
timescale for this, given the clump masses and total 
masses of our models, is typically around 500 Myr (see 
Bournaud et al. 2007). This is slower than gaseous tur- 
bulence dissipation, and therefore not as important for 
overall shrinkage when the gas fraction is large. 

Contraction also occurs if we consider the size of a disk 
to be determined by angular momentum. Baryons in the 
central few kpcs of interacting galaxies loose significant 
angular momentum, while baryons initially in the outer 
disks carry away a large fraction of t he initial angula r 
momentum in tidal tails (see review in IBournaudl[20T0l) . 
At first order, the angular momentum of baryons near 
the corotation region, which is typically about the half 
mass radius (a few kpc for the disk masses considered 
here) is be roughly conserved. 
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For a mass AI encompassed by a radius i?, and a cir- 
cular velocity V = \J GM / R, angular momentum con- 
servation implies that: 

i?cxl/FMoc 

If gas initially outside the half-mass radius dissipates 
its turbulent energy and moves towards the center, the 
central mass increases, possibly as much as a factor of 
2. Then the half-mass radius can decrease by a factor of 
8, all the while conserving angular momentum. Further 
shrinkage can be caused by angular momentum removal 
through tidally expelled baryons. Only this tidal removal 
is present in models with a smooth, non-turbulent ISM. 

3.4. Star formation history 

The star formation histories of mergers CI, GIF and 
SI are shown in Figure [51 The efhciency was cali- 
brated in each model to result in a SFR of approximately 
100 Mq yr~^ in each pre-merger galaxy, typical for z ^ 2 
disk galaxies with similar masses. The star formation 
rate during the merger reaches a factor of '-^lO higher 
than in the pre-merger pair in all cases, with peak val- 
ues of 1000-2000 Mq yr~^ This rate is consistent with 
the SFRs estimated in S ubMillimeter Galaxies (SMGs; 
e.g.. iTacconi et alll2008D . The star formation peak can 
occur earlier-on or later-on in various models, but the 
maximal SFR and the duration of the starburst remain 
comparable, leading to globally similar amounts of gas- 
to-star conversion for the cooling and stabilized models. 
The fraction of remaining gas at the first pericenter and 
300 Myr later are given in Table 2: no major or system- 
atical differences are found between the cooling and the 
stabilized ISM models for any given orbit. 

Both the cooling and the stabilized ISM models reach 
rates of star formation that are realistic for observed 
z ^ 2 disks and starbursting mergers. The rates are 
also consistent with each other. This does not mean 
that gas cooling below 10^ K in a clumpy turbulent 
ISM has no effect on star formation in mergers. High- 
resolution models explicitly including gas cooling could, 
in principle, resol ve the densest gas ph ases that are actu- 
allv star-forming (iTasker fc Tanll2009 l: IAgertz et al.ll2009l: 



merger disk pair (with both the cooling and stabilized 
EoS). This factor is somewhat higher than the typi- 
cal value in larg er samples of merger simulations (e.g . 
Pi Matteo eraIT[20 07. 2008; iMartig fc BournaudI mm 
Stewart et al ]l008) but not in large disagreement. This 
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iBournaud et al.l l2010[) . but models using an artificially 
thermalized ISM do not reach such high gas densities and 
they have to describe star formation as a low-efficiency 
process pervasively distributed throughout a low-density 
gas. The calibration of the efficiency is subjective when 
dense clouds are not made self-consistently by the model. 
Here we choose an efficiency in the stabilized ISM model 
that gives the same star formation rate as the cooling 
model for an isolated disk. Had we used the same pre- 
scription for star formation in each case, then the cool- 
ing models would have had much higher star formation 
rates and their collisions would have been more gas-free. 
While we compare the cooling and stabilized models with 
similar SFRs and resulting gas fractions, the EoS used to 
model cold gas imposes most of this gas to be relatively 
cold < 10'' K. Different feedback models might result in 
higher amounts of hot extended gas, which could poten- 
tially cool down later-on and accrete onto a re-formed 
disk (see discusion in Sections 3.4 and 3.5. 

Globally, the star formation rate in our merger mod- 
els peaks at about ~ 10 times the value of the pre- 



factor of about 10 also appears somewhat higher t han in 
observations of merge r s and close pairs (Bergvall et al 



ODseryations oi merge r s ana close pairs t ijergvaii et al 
2001 iBell et all [20051: iJogee et al.l 120091: iRobaina et al 



20091) . However, these observations indicate the SFR en- 
hancement at a random instant of an interaction, and 
what they consider to be major mergers are not strictly 
equal-mass ones. The factor 10 at the peak SFR in 
our models, which is reached only for a short period, 
is broadly consistent with the factor of 3-4 indicated by 
these observational studies. We did not here focus on 
particularly efficient cases where the SFR enhancement 
in a merger can reach factors larg er than a hundred (as in 
,Mihos &: Hernauistl996. or .Springel fc Her n auist.2005.) , 
but rather on random, representative orbits. 

The peak in the merger-induced star formation rate 
begins somewhat earlier in the cooling models than in 
the stabilized ones, just after the first pericenter passage 
(Fig. |6]) . This earlier gas consumption cannot cause the 
lower disk fraction of the merger remnants because (i) ii 
occurs only once the interaction process starts, so the gas 
fractions during the interaction remain quite similar in 
both models (Table 2) and (ii) the increased feedback in 
model GIF delays the peak of star formation until a bit 
later than in model SI, but the final merger remnant is 
almost as compact and disk-free as for the initial model 
Gl. 

The earlier merger-induced star formation activity in 
the coolin g models is c onsistent with the results ob- 
tained by iTevssier et al. (2010) in low-redshift merger 
models (see also iSait oh et al.l l2009l for resolved star- 
forming regions in early pahses of galaxy interactions). 
High-resolution models with a low temperature fioor can 
capture ISM turbulence, local dense shocks, and frag- 
mentation into star-forming clouds. Stars form in these 
clouds when they are still in the main disk, long be- 
fore they have accreted to the galaxy center. Such mod- 
elling is required to reproduce correctly the extent of star 
formation in early- and late-stage interactions, which 
is ni i ssed by standard pre scriptions (see also iBarnesI 
120041 iGhien fc Barnes! I2010t ). On the other hand, low- 
resolution or stabilized ISM models get intense star for- 
mation only in the dense centralized gas that follows the 
merger-driven inflow. Dense clouds do not form in the 
disk. The delay in star formation for this stabilized case 
is comparable to the rotation period of the entire galax- 
ies. 

3.5. Feedback and clump evolution in mergers 

Our simulations include a kinetic model for the en- 
ergy feedback from Type II supernovae. Other sorts of 
feedback such as Type la supernovae, stellar winds, and 
mass return are important for high-rcdshift galaxy evolu- 
tion ( Agertz et al. 2010: Martig & Bournaud 2010), but 
the timescales for this feedback are typically longer than 
the timescale for a major merger of massive disk galax- 
ies. Another sort of feedback with a short timescale 
is radiation pressure from young massive stars. Gi- 
ant clumps in high-redshift, gas-rich galaxies could be 
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Fig. 6. — Star formation history for the merger models along orbits 1, 2 and 3, with various EoS and feedback parameter. The pericenter 
passage corresponds to i = Myr in each case. 



sufficiently massive to self-regulate their star forma- 
tion so that they "survive" feedback for several hun- 
dred Myr ([Krtimholz fc Dekel 20 .10). This however re- 
mains de bated (iMurrav at al".ll2'o'icjl) . The "survival" defi- 
nition of iKrumholz fc Dekell ()2010l ) implies that a bound 
stellar clump remains but not necessarily that the gas 
is retained by this clump. Expelled gas could be re- 
accreted by the clump and start forming stars again 
(|Pfiamm-Altenburg fc KrouD a''2009^. 

These additional feedback mechanisms are not present 
in our models. Nevertheless, we found that varying the 
efficiency of the supernovae feedback (in runs GIF and 
LM-C2F) did not result in major changes to the prop- 
erties of the final ETGs. Models CI and GIF also have 
similar line-of-sight velocity dispersions in their gas com- 
ponent (Fig. [3] and S]). The main source of turbulence in 
high-redshift disk galaxies is generally considered to be 
gravitational energy released thro ugh clumping instabili- 
ties a n d/or inward mass a ccretion (jElmegreen fc BurkertI 
120101 iDekel et all l2009| ). Observations have sug- 
gested that the local velocity disperion scales with the 
local gas density and/q r star fo rmation surface d en- 
sity (jLehnert et al.1 120091 see also lGreen et al.l 120101 for 
low-redshift analogues) which could suggest energy in- 
put from stellar feedback. The local velocity disper- 
sion scales with the local surface density also in models 
dominated by gravity -driven turbulence such as those in 
iBournaud et al.|[2009l Gas turbul ence can be further in- 
creas ed by tidal forces in me rgers (lElmegreen et al.| | 1993 : 
Irwinlll994HElinegreen et al.lll995l: IBournaud et al-lboOS : 
Teyssier et al.l 120101 ). Stellar feedback is thus not ex- 
pected to be the main energy source of the turbulent 
motions. The sa me ap pears to be true in local galaxies 
(jBournaud et al.|[2010[ ). 

While the pre-merger disk galaxies are already clumpy, 
the main gas clumps during the merger tend to be denser 
and more massive (see for instance the time sequence 
on Figure 1): some are pre-existing clumps that accrete 
more mass, some new clumps form with high masses al- 
lowed by the high gas densities and velocity dispersions. 
This maintains a clumpy stellar morphology during and 
after the merger. Most clumps are nevertheless not long- 
lived throughout the merger process. The average stellar 
age of the five largest stellar clumps in the final ETG in 
model Gl was measured 450 Myr after the pericenter pas- 
sage, which is about 200 Myr after the final coalescence. 
These ages range from 80 Myr to 175 Myr. Hence these 
clumps are relatively young compared to the merger and 
formed in the latest phases. This suggests that our re- 
sults do not directly rely on a low feedback efficiency, but 



primarily result from a strongly turbulent ISM, in which 
massive clumps naturally arise. 

We have not explored the modeling of other sources 
of stellar feedback (radiation pressure, winds, etc). We 
already discussed models with an efficiency up to 100% 
in our supernovae feedback scheme, but Agertz et al. 
(2010) suggested that higher efficiencies, up to 500%, 
could be more realistic - presumably to account for other 
sorts of feedback that are not modeled, such as pho- 
toionization, or radiative feedback, the role of which at 
high redshift remains debated (e.g. [Murray et al.l 120101 : 
IKrumholz fc Dekell [20Tol: IGenel et aLn2010l) . Thus, we 
performed tests with SN feedback efficiency increased 
to 500%, for the cooling and stabilized models on or- 
bit 1. The results (see Table 2) show some variations, 
as stronger feedback fuels more extended and hotter gas 
reservoirs in the early phases of the merger, which adds 
mass to a re-formed disk component in the post-merger 
phases (following the process discussed, e.g., by Springel 
& Hcrnquist 2005 and Governato et al. 2009). Nev- 
ertheless these variations are relatively minor compared 
to the initial discrepancies between models Gl and SI, 
both in terms of final disk fraction and final compact- 
ness. Furthermore, the differences between models GIF 
and SIF, or G1F5 and G1F5 respectively, are about of 
the same amplitude as the differences between the initial 
models Gl and SI. This suggests that, while a stronger 
feedback can moderately increase the disk fraction and 
half-mass radius of the final ETG, this effect is decoupled 
from the impact of using a cooling model generating tur- 
bulent gas rather than a thermally stabilized model. It 
does not strongly impact our previous conclusions on the 
effect of resolving gas cooling and strong ISM turbulence 
on the stellar size evolution and disk fraction in the final 
merger remnants. 

3.6. Mergers of lower-mass galaxies 

Our main set of simulations used galaxy models that 
were relatively massive (without being unrealistic) for 
the typical disks and spheroids observed at z ~ 2. 

Lower- mass disk galaxies at z ~ 2, with circular ve- 
locities of 100-200 km s~^, also show giant clumps of 
star formati on and gas turbulent speeds of several tens 
of km s-i (jForster Schreiber et all 120091 : iWright et ahl 
l2009t) . like the pre-merger galaxies in our models. ISM 
turbulence and dumpiness should thus have qualitatively 
similar effects in high and low mass major mergers. This 
can be explicitly checked using models LM-G2, LM-S2 
and LM-G2F. 

The structural parameters of these low-mass mergers 
remnant are given in Table 2, and the edge-on stellar 
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Fig. 7. — Outcome of the lower-mass merger models LM-S2 and 
LM-C2 after full relaxation of the merger remnant. The projected 
stellar surface density is shown in each case, with an edge-on ori- 
entation with respect to the stellar spin axis. Although a weak 
disky distortion is observed in the cooling model LM-C2, only the 
stabilized ISM model LM-S2 has a massive and extended rotating 
disk component, and it also has a larger half-mass radius. 

distributions shown in Figure [71 The resuhs are simi- 
lar to what has been detailed previously for mergers of 
higher-mass galaxies: the gas velocity dispersions in the 
merging phase are higher in the cooling model, and the 
merger remnant is significantly more compact, with only 
a low fraction of its baryons in a rotating disk, while 
a significant disk component is found in the stabilized 
ISM model. In detail, the relaxed ETG after the merger 
is comparatively more compact, with a stellar half-mass 
radius reduced by a factor of 2.5 in the cooling models 
compared to the stabilized ISM model, versus 2-2.2 for 
the more massive cases. 

4. COMPARISON TO OBSERVED GALAXY POPULATIONS 

4.1. Kinematical signatures and dispersion- dominated 

systems 

The mergers modeled in realistic high-redshift con- 
ditions develop very high gas velocity dispersions, and 
could thus resemble some dispersion-dominated systems 
(DDS) observed in z ^ 2 surveys. To further probe 
whether such mergers would be easily identified in obser- 
vations, three of us experienced in high-redshift spectro- 
scopic data (TC, BE, KS) performed an eye classification 
of the velocity and dispersion fields from Figure 3 with 
gas density maps, and from Figure 4 without the den- 
sity maps. The test was done without knowledge of the 
nature of the modeled systems. This test aims at under- 
standing whether the kinematics of the simulated merg- 
ers would be recognized as typical for major mergers by 
observers familiar with the typical signatures expected 
for mergers. Additional instrumental/observational ef- 
fects such as limited resolution, sensitivity, and noise, 
were not accounted for: hence this test relates only to 
the physical properties shown on Figures 3 and 4. Real 
observations could in addition miss small-scale or low- 
fiux signatures and make the identification of mergers 
even more uncertain. 

Model I is always recognized as an isolated disturbed 
disk, with some suggestions that the top-left clump could 
be a minor merger. The stabilized-ISM merger models 
are most often identified as mergers (two-thirds of the 
votes), generally with a high confidence level that these 
are binary major mergers. The cooling merger models 
were classified as DDS of unclear nature (about one-third 
of the votes), DDS of merger origin (one-third), or merger 
(one-third). More than half of the votes suggesting a 
merger or DDS merger also suggested a group or multiple 
merger rather than a binary major merger. This confirms 
that (i) on-going mergers with clumpy turbulent gas 



could be classified as DDS, and (ii) typical high-redshift 
mergers have substantially different kinematics from low- 
redshift mergers or models with artificially stabilized gas, 
which have smoother velocity gradients along extended 
tidal tails. 

Many observed DDS have small sizes, so the high ob- 
served dispersions could sometimes result fro m blending 
at lo w resolution, but seem real in general ()Law et al.l 
l200a . There are nevertheless DDS that are relatively 
massive, extended, with chaotic velocity fields. Exam- 
ples incl ude Q1623-BX453j QITOOO-BXTIO, and others in 
ILaw et all (|2061 and in iF6rster Schreiber et al.l (|2009D . 
These could be major mergers of initially clumpy and 
turbulent galaxies, but their merger nature has not been 
clearly recognized yet because they lack the usual known 
signatures of mergers. 

Disk galaxies with lower masses tend to be even more 
tur bulent compared to the rotat ion speed (lower V/a, 
e.g. iForster Schreiber et aLll2009D . Reaching dispersion- 
dominated phases in mergers of such galaxies should be 
easier at lower mass. This might explain why observed 
DDS generally have small stellar masses, even though 
they are too numerous to all be major mergers. Over- 
all, DDS are often sugges ted to result from violent gassy 
collapse ([Law "eFall [20091 ) ■ which our models show can 
occur in high-redshift mergers. 

The local turbulent speed of gas in our cooling mod- 
els is typically 150-200 km s~^ (Table 1), but the pro- 
jected line-of-sight dispersions can be up to 500 km s~^ 
for central clumps and up to 350 km s^^ for other clumps. 
These very high dispersions arise on short phases, shortly 
after the interaction pericenter and about at the peak of 
SFR. These dispersions may seem higher than what is 
found in high-redshift spectroscopic surveys. Note how- 
ever that our models CI to C3 involve galaxies that 
are more massive than most objects targeted in spec- 
troscopic surveys, so they naturally have higher circu- 
lar, and higher dispersion velocities at fixed V/a. Fur- 
thermore, these high dispersions would be blended out 
by seeing-limited observations without Adaptive Optics. 
Substantial samples of Adaptive Optics spectroscopic 
data really exist only for disk galaxies (e.g.. lGenzel et aLl 
1201 Ih . Other high- resolution data are unlikely to have 
targeted nearly equal-mass mergers at the peak of their 
high velocity dispersion phase. Whether or not the peak 
gas dispersions in our models are too high for real 2; ~ 2 
objects is thus presently unknown; if this is the case, 
the effects of turbulence dissipation and size reduction in 
major mergers might be somewhat more modest. 

4.2. Submillimeter galaxies as on-going mergers? 

SMGs are among the most actively star-forming ob- 
jects at high redshift, and hence generally proposed to 
be t he most actively sta rbursting phases of major merg- 
ers (jTacconi et al.l [20081 ). An interesting test is then to 
compare the gas dynamics in our models to observations 
of SMGs. The properties of a representative z = 2.3 
SMG, SMMJ2135-0102, have been map ped spatially ow- 
ing to gravitational magnification by ISwinbank et alj 
(|2010f ). While the global properties are reminiscent of 
a gas-rich major merger, star formation takes place in 
a few giant clumps with sizes of 100-10 00 pc, as in our 
merger models (Fig. [ij. Interestingly, ISwinbank et al.l 
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Fig. 8. — Integrated spectrum of model CI, 55 Myr after the 
pericenter passage, i.e., the second snapshot on Fig. 1, within a 
central aperture of 4 kpc. Gas denser than 500 cm~^ was selected 
as a proxy for dense CO-emitting molecular gas. The dashed spec- 
trum is for lower density gas (denser than 10 cm"-^). The multi- 
component aspect of the dense gas spectrum is reminiscent of the 
CO spectra observed in SMGs (see text for details). 

()2010[) and lDanielson et all (|2010f ) noted that the molec- 
ular gas spectra have a multi-component shape. The 
spectrum is more complex than the usual double-horn 
profile observed for rotating gas disks (including clumpy 
high-redshift disks, Daddi et al. 2010), and is best 
fitted by a combination of four main emitting com- 
ponents. Multi-component spectra in SMGs that are 
not just doub l e-horn profiles were also suspected by 
iKnudsen et"all (|2009fl . 

We modeled a "molecular gas" spectrum for model CI, 
55 Myr after the pericenter passage, using gas denser 
than 500 cm~^ as a proxy for molecular gas. The mod- 
eled spectrum, corresponding to the second snapshot on 
Figure 1, is shown in Figure |8l It has a multi-component 
aspect similar to that in the observation. Our models did 
not aim to reproduce this SMG in particular, and each 
individual component is not reproduced, but the ampli- 
tude of the e mitting peaks is consistent with the spectra 
observed by iPani elson et al.l (|2010[ ) . This property re- 
sults from the dense cold gas (and star formation) being 
mostly located in a few giant clumps. Interestingly, the 
stellar mass at this instant in model CI is 2 x 10^*^ M0, 
with a gas fraction of 40%, close to the estimates for 
SMMJ2135-0102. The star formation rate in the model 
is 1400 M0 yr^^, some what higher than o bservational es- 
timates for this object (jivison et al.ll2010l ) but not incon- 
sistent with SMG-like activity in general. The modeled 
molecular gas spectrum has a full width at half maximum 
and a full width at zero intensity that are quite close to 
the observed values, suggesting quantitative agreement 
on global gas dynamics between our model and a typi- 
cal SM G of similar mass. Furthermore Danielso n et al.l 
(j2010( ) noted that the denser molecular gas tracers (such 
as high-J CO lines) have more peaked spectra, while 
lower-density tracers have smoother spectra (Greve et 
al. 2005 also found relatively smooth low- J spectra in 
SMGs). This is consistent with the peaked components 
being associated with dense star-forming clumps, and 
with our models showing smoother spectra for lower- 
density gas (Fig. El). Higher- resolution simulations re- 
solving cooling down to lower temperatures would be re- 
quired to qualitatively explore further this temperature 
dependence of molecular gas dynamics in SMGs. 

Standard models of high-redshift mergers using 



warmed stabilized gas to model the ISM are only 
mildly successful in explaining the properties of SMGs. 
They can model 805/im fluxes that account for SMGs 
(jNaravanan et all 1201 0). but only when the pre-merger 
disks are calibrated to high submillimeter iluxes: the 
typical enhancement of the 850/im flux during the 
merger/SMG phase is not larger than a factor of 2- 
3. These models also have smooth gas kinematics. If 
dust follows the scattering of gas in a three-dimensional 
clumpy distribution in and around the centrally collapsed 
system, as in our models, then the obscuration of star for- 
mation co uld be very high unde r some lines of sight, like 
in SMGs (|Thronson et al.l[T990l) . 

Another model of SMGs suggests that they are the 
high-m ass end of cold-flow galaxies, with no mergers in- 
volved ()Dave et al.l 120101 ). Some SMGs could be like 
this, particularly those with large radii, relatively low 
star-formation rates, and relatively long durations com- 
pared to the orbit time. The SMGs in our model have 
much stronger bursts and the disk gas is consumed more 
quickly than in the cold flow models, so we expect a 
different proportion of elemental abundances in these 
two cases, and different star formation histories. Multi- 
component gas spectra could be a confirmation of a 
clumpy merger origin for SMGs, if they appear to be 
commonly observed in larger samples. 

4.3. Clumpy morphologies in high-redshift mergers and 

ETGs 

The stellar distribution in our models shortly after the 
pericenter resembles a single "clumpy galaxy" , but it is 
more asymmetric than the pre-merger clumpy disk mod- 
els (see second and third panels in Figure 1). This could 
be consistent with t h e pop ulation of "assembly galaxies" 
in lElmegreen et al.l (|2007f) . with entirely clumpy mor- 
phologies not consistent with their being isolated disks. 
Other e xamples of high-re dshift clumpy mergers can be 
found in lLotz et all (|2006[) . 

Our modeled high-redshift clumpy mergers lack very 
prominent tidal tails. There are long tidal features, but 
they are less contrasted than the star-forming clumps. 
This differs from low-redshift mergers, where long tidal 
tails wit h ample and smooth vel ocity gradients are typi- 
cal fe.g.. lBournaud et al.ll2004bl ). 

Once the mergers are relaxed, stellar clumps are seen 
in the central kpcs of the final ETGs (Fig. ^. They 
are relatively young (1 — 2 x 10® yr in the final snap- 
shots) and their masses are all lower than 2 x 10® Mq 
- while our pre-merger disk models harbor clumps of 
several 10® M0. A lower clump mass is expected from 
stabilizati on of the merger r e mnant by the ma in stellar 
spheroid f Martig et al.l [20091: iDekel et al.ll2009[ ). Resid- 
ual star formation is associated with those clumps after 
the merger, consistent with high-redshift red, but not 
completely dead, ETGs (jTonini et al.ll2010l ). 

These clumps might be harder to observe than the 
most massive ones in z ~ 2 star-forming disk galaxies, 
but young stellar clumps have actually been observed in 
the central body of high-r edshift ellipticals in the Hubble 
Ultra Deep Field (UDF; lElmegreen et al.ll2005l see ex- 
amples in Fig.|n|). Those clumps have young stellar ages, 
and their masses are typically a fraction 10~'^-10~'* of 
their host galaxy masses - lower than in clumpy disks but 
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Fig. 9. — Four examples of elliptical galaxies iu the UDF where 
unsharp masking unveils clumps in/around the main body of 
the ellipticals. Other examples and details can be found in 
(|Elmegrecn et al. 2005). The clumps are relatively young, and less 
massive than the clumps in UDF disk galaxies, consistent with the 
clumps in early- type galaxies in the outcome of our merger models. 

consistent with our ETG models. Other clumps are also 
seen around several ETGs in the Mancini et al. (2010) 
sample. 

The outermost clumps, those formed at more than 
50 kpc from the center in model CI, could become dwarf 
galaxies formed diiring t he merger (tidal dwarf galaxies 
(TDGs) iDuc et a l.|[2000t). This would lend support to 
recent suggestions that the rate of TDG form ation could 
be high at high redshift (|Kroupa et al.l[2010l ). 

4.4. Compact high-redshift ETGs 

Our merger models with cooling predict final stellar 
distributions that are relatively compact. The proper- 
ties of the final ETGs in our cooling models are com- 
parable with the properties of the three most co mpact 
high-redshift ETGs in the iMancini et"an ()2010l ) sam- 
ple with similar masses, as well as some objects in 
iCappellari et al.l (|2009D . Faint extended stellar halos 
with high Sersic indices are present in the observations, 
like the one surrounding the massive kpc-sized core in 
the final ETG of model CI (Fig. 4). 

Real mergers should involve various mass fractions, 
gas fractions, and degrees of stability and dumpi- 
ness in disk galaxies. Thus, they should form a va- 
riety of ETGs, ranging fro m very compact objects 
(a s spectroscopically confirmed. Ivan Dokkum et al|[200l 
ICappellari et al. 2009) to more extended ones (as present 
in the Mancini et al. sample). If high-redshift galax- 
ies with lower masses are relatively more turbulent 
(jForster Schreiber et al.ll2009l ). then their mergers could 
be even more dissipative than the more massive disk 
galaxies. Then the most compact ETGs would be the 
lowest-mass (Section 3.5). Such a trend is observed 
among the various samp les of high-redshift spheroids (see 
F ig. 4 inlM ancini et al. 2010( ). 

iWuvts et al., (,2010,1 also proposed recently that high- 
redshift mergers lead to the formation of compact ETGs. 
Their simulations have a stabilized, pressurized ISM 
model (although implementation details differ from our 



"stabilized" models). There is some compacting of the 
system during the merger, but at the same level as 
in our stabilized models and not at the level of our 
cooling models of clumpy galaxy mergers. Actually, 
to fo rm ETGs th a t are about as compact as 2; ~ 2 
ones, iWuyts et al.l ()2010l ) had to start with pre-merger 
disk galaxies that were twice as compact as our ini- 
tial disk galaxies, with a disk scale length half as large. 
Which initial conditions are more realistic in terms of 
disk sizes is unknown. The known size evolution be- 
tween z ^ 2 an d z — (Fr anx et al. 2008) mostly con- 
cerns spheroids. iBouche et al.l (|2007l ) propose that high- 
redshift disks have the same velocity-size relation as disk 
galaxies at z = and a spin parameter 2 to 6 times larger 
than in the Wuyts et a l. initial cond i tions, but this was 
called into question bv lDutton et al.l ()2010t ). Our initial 
disk galaxies still seem realistic for massive star-forming 
disks at z ~ 2: a radial scale- length of 5 kpc for a circular 
velocity of 245 km s~^, to be c ompared for instance to 
the su b-sample of disk galaxies in lForster Schreiber et al.1 
(|200a Fig.18). 

It thus appears that a major merger model can pro- 
duce the transformation of a typical z ^ 2 star-forming 
disk galaxies (hence relatively extended) into a typical 
high-redshift ETC, only if ISM cooling, turbulence and 
dumpiness are captured in the simulation. Observations 
do suggest that compact ETGs formed from relatively 
diffuse disk s, rather than disk ga laxies that are already 
compact (Ricci ardelli et al"1l2010t ). It nevertheless seems 
that our models predict Sersic indices that are somewhat 
higher than those estim ated for compact h igh-redshift 
ETGs. In observations, iPuech et al.l ()2006[ ) also noted 
that compact luminous galaxies at high redshift rarely 
show rotational support, but instead have large disper- 
sions and chaotic velocity fields. They suggested that 
high-redshift mergers produce compact remnants. 

4.5. Long-term evolution towards low-redshift ellipticals 

High-redshift ETGs probably do not survive in the 
form of compact objects in the nearby Universe. They 
should grow mostly throug h minor rnerger s, which can 
increase their stellar sizes ( Naab et al.ll2009l) . In this dry 
process, the small and compact disk components remain- 
ing after high-redshift mergers (Fig. 1 at t = 625) could 
undergo inefficient star formation (|Martig et al.l l2009| ) 
and persist as compact disks of residual gas and stars, 
as often found in the central kpc of n earby ellipticals 
(|Young et al.ll2008HCrocker et al.ll20Tol) . 

Of course some of these high-redshift ETGs could ac- 
crete external gas in sufficient amount to reform a bright 
and massive stellar disk, in a "spiral rebuilding " scenario, 
as proposed for instance bv lPuech et al.l ()2006t ). However 
our models show that such disk rebuilding does not hap- 
pen spontaneously from gas expelled during the merger 
(in tidal tails and/or through feedback processes) that 
would later on fall back onto a disk. Large accretion of 
external gas would be required to rebuild a massive star- 
forming disk galaxy. Models bv lMartig et alj (|2009[ ) have 
shown that this can naturally happen in some haloes in 
A-CDM cosmology, when the late mass assembly history 
is dominated by smooth infall rather than mergers. 

5. CONCLUSIONS 
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The ISM is very turbulent and clumpy in high-redshift 
galaxies. Simulations of high-redshift mergers with cool- 
ing and sufficient resolution to develop turbulence and 
dense gas clouds show major differences compared to 
models with smooth, thermally supported disks: 

• On-going mergers of representative high-redshift 
galaxies can have very irregular gaseous spheroids 
during the merger, with clumps and filaments 
throughout, instead of rotating gas disks sur- 
rounded by long tidal tails with ample and smooth 
velocity gradients, as is common in low-redshift 
mergers. 

• The gas velocity dispersion in these mergers be- 
comes very high. At some stage, most of the gas 
kinetic energy is in the form of a three-dimensional 
velocity dispersion. Such mergers could take the 
appearance of dispersion-dominated galaxies with 
high gas dispersions and chaotic velocity fields. 
Their spectral properties are also consistent with 
SMGs. 

• Turbulent and clumpy gas undergoes a violent and 
dissipative collapse in mergers. The final ETGs 
in these mergers are relatively compact, with high 
Sersic indices and faint outer stellar halos. 

• The masses and sizes of disk components that sur- 
vive or re-form after major mergers are strongly 
reduced when the hydrodynamics of the cold tur- 
bulent ISM is taken into account. Large disks at 
high redshift should form mainly through cold ac- 
cretion rather than mergers. 

• Strong feedback processes can somewhat increase 
the size and disk fraction of the final ETGs formed 
in these mergers, but the impact of directly model- 
ing a cold turbulent ISM phase rather than using a 
thermally pressurized model remains independent 
from the efhciency of stellar feedback. 

Wet mergers at high redshift should involve galaxies 
that have turbulent and clumpy gas-rich disks. Merg- 
ers of such galaxies can undergo dispersion-dominated 
phases. They can also produce starbursts with SMG- 
like properties. The late phases are accompanied by dis- 
sipative collapse and the formation of compact ETGs. 



This is consist ent with the evolutiona ry sequence pro- 
posed, e.g., bv iRicciardelli et all ()2010[ ). Massive rotat- 
ing disk components cannot easily survive such mergers, 
a property that is not reproduced correctly by simula- 
tions where ISM turbulence and dumpiness are unre- 
solved and modelled with a thermally pressurized equa- 
tion of state. Large rotating g as-rich disks a t high red- 
shift, such as those observed bv lGenzel et al.l ()2006D and 
others, are probably not the outcome of a mass assem- 
bly dominated by mergers. Large far-outer reser- 
voirs of gas, pre-existing to mergers or fed by strong 
feedback during the mergers, coul d participate to some 
disk rebuilding (as observed in the lGovernato et al.ll2009l 
model). Stellar mass- loss and gravitational torquing 
over cosmological timcscalcs can also help preserve large 
disks (as in theLM artig & Bournaud 2010 model, see also 
iLeitner fc Kravtsovi 2010), but continuous cosmological 
infall of fresh gas probably dominates the assembly of 
disk-dominated galaxies down to z — 0. 

If gas-rich mergers at high-redshift are characterized by 
strong turbulence, the spread of the gas density distribu- 
tion could be very large, as it scales with the turbulent 
mach number (Krumholz & Thompson 2007). Not only 
the average gas density would be high, but there would 
also be an excess of dense gas at number densities of thou- 
sands per cm^ and more, compared to non-interacting 
systems. Dense gas excess has alre ady been found frorn . 
low-redshift merger models (e.g., lJuneau et al.l 120091 : 
iTevssier et al.l 120101 ) and could explain the enhanced 
HCN/CO ratios in nearby ULIRGs. An excessive frac- 
tion of dense gas could also increase the star formation 
rate at fixed average gas density, possibly consistent with 
recent observations of star formation in disks and merg- 
ers suggesting enhanced efficiency of gas consumption in 
mergers (Daddi et al.ii2010k .Genzel et al.ii201Q ). 
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